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Abstract: 

An exciton gas on a lattice is analyzed in terms of a convergent hopping expansion. For a given chemical 
potential our calculation provides a sufficient condition for the hopping rate to obtain an exponential 
decay of the exciton correlation function. This result indicates the existence of a Mott phase in which 
strong fluctuations destroy the long range correlations in the exciton gas at any temperature, either by 
thermal or by quantum fluctuations. 


1 Introduction 

Coupled quantum wells represents a class of systems which allows us to study strongly interacting particles 
under controllable conditions. They are conceptually simple: negative electrons are trapped in a two- 
dimensional plane, while an equal number of positive holes is trapped in a parallel plane a distance D 
away (see Fig. H}. One of the appeals of such systems is that the electron and hole wavefunctions have 
very little overlap, so that the excitons can have very long lifetime (> 100 ns), and therefore, they can 
be treated as metastable particles to which quasi-equilibrium statistics applies. 

The exciton gas is effectively a hardcore Bose gas, whose groundstate should be a superfluid at low 
density or a Mott state at high density. The latter requires a lattice such that a lattice commensurate 
state can be formed. A superfluid state of excitons in coupled quantum wells was predicted some time 
ago in Ref. pf. Several subsequent theoretical studies [2H93 have suggested that superfluidity should be 
manifested as persistent electric currents, quasi-Josephson phenomena and unusual properties in strong 
magnetic fields. In the past ten years, a number of experimental studies have focused on the observation of 
the superfluid behavior (101 - 1171 . The transition from an exciton gas to an electron plasma in GaAsGaAlAs 
quantum wells was analyzed in the framework of many body effects, considering the dynamical screening 
of the Coulomb interaction in the one-particle properties of the carriers and in the two-particle properties 
of electron-hole pairs [15] . This was also studied for a 2D electron-hole system, considering the exciton 
self-stabilization mechanism, caused by the screening suppression due to the exciton formation [19] . We 
recently studied a double-layer exciton gas in mean field approximation and found a transition to a Mott 
state at high densities [20] ■ Although the destruction of the superfluid state by strong collisions in the 
dense exciton gas is plausible, a mean field approximation is not a very trusted tool to analyze a two- 
dimensional system. Therefore, it remains to be shown by an independent and reliable method beyond a 
mean field approximation that the long-range superfluid correlations are destroyed in a dense system for 
all temperatures. 

Since an exciton is a strongly bound pair of an electron and a hole, we assume that these pairs cannot 
dissociate or transform into a photon by recombination. This can be justified by a strong Coulomb 
interaction and a sufficiently short time scale that is shorter than the recombination rate. Moreover, 
we implement a lattice structure in the layers to allow the excitons to form a commensurate Mott state 
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Figure 1: Two coupled quantum wells with lattice gates, realized for electrons (filled circles) in one well 
and holes (empty circles) in the other. 


by filling each well of the lattice with an exciton. A lattice can be realized by an electrically charged 
gate which is periodically structured [5T] . The gate can also be used to control the density of excitons 
via a chemical potential /i. Then a pure exciton gas in the periodic potential can be described by the 
Hamiltonian 

H = — + /-iS rr ')ala r ' + h.c. (1) 

rr' 

where the sites r, r' are the minima of the potential wells and J rr > is a nearest neighbor hopping rate: 


f J if r, r' are nearest neighbors 
10 otherwise 


The lattice structure is characterized by the number of nearest neighbors c (connectivity). The exciton 
creation operator aj is composed of the electron creation operator c\ r and the hole creation operator c' h 
as 

al = r . ( 2 ) 

The form of the exciton operator in Eq. © implies that the excitons can tunnel freely with the hopping 
rate J rr > under the restriction that they are composite particles which obey the Pauli principle. In other 
words, at most one exciton can occupy a site of the lattice. Therefore, it is a hardcore boson, similar to 
the exciton described by the effective Hamiltonian of Ref. [55]. 


2 Conditions for a Mott phase 


Starting from the Hamiltonian © we consider a grand canonical ensemble of excitons with the partition 
function Z = Tre~^ H at temperature T (/3 = l/ksT). The trace Tr is taken with respect to all exciton 
states. The exciton correlation function then reads 

(ef* tH ale-l> tH ef* H a r .e-l* H ) = ^-Tre^ 13 e^ tH e^' H a r ,e~^ H (3) 

Zj 

in imaginary time representation with 0 < t — t' < 1 [23 . 

The partition function Z of a one-dimensional hard-core Bose gas is given as the determinant of 
non-interacting fermions (cf. Ref. [24] ). The reason for the equivalence of non-interacting fermions and 
hardcore bosons in one dimension is the fact that fermions cannot exchange their positions but must 
obey the repulsive Pauli principle. In contrast to a non-interacting Bose gas, this system exhibits two 
Mott and one an intermediate incommensurate phase. Using this example we can expand its free energy 
F = —/3 _1 log Z in terms of the hopping matrix J as 


F = log(l + e^e 133 ) = -TrJ - ilY[log(l + e^*)] 


R / 


;>i 


1 


1 + 


(e-^ J - 1) 


(4) 
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Figure 2: Phase boundary for the dimensionless chemical potential fl/j, vs. the dimensionless hopping 
flj of excitons on a square lattice: A Mott phase with exponentially decaying exciton correlations exists 
above the solid curve. The dashed line is the mean field condition flj < flp, for a Mott phase. 


provided that J is small in comparison to /i. For flj 1 (strong quantum fluctuations) the inequality is 
valid for p, < J. This example gives us an idea about the convergency of the hopping expansion in higher 
dimensions. 

For a dilute two-dimensional system we expect a formation of a superfluid state with power law 
correlations {e^ tH a\e~^ tH e^' H ar'e - ^' H ) ~ \r — r'\~ a [25]. On the other hand, for sufficiently high 
density the long-range phase correlations will be destroyed by frequent collisions of the excitons. This 
effect leads eventually to an exponentially decaying correlation function leP tH a\e~^ tH e ^ 4 ff a r /e _/3 * H ) ~ 
exp(— \r — r ’|/£) with the correlation length £. Such a state is either a thermal exciton gas with a fluctu¬ 
ating density at high temperatures or a Mott state with fluctuating phases but non-fluctuating density 
at low temperatures. These two states change from one to the other upon reducing the temperature. If 
the two phases are connected by a crossover or by a phase transition is not clear at this point. However, 
both phases are characterized by a gapped excitation spectrum. A mean held approximation indicates 
a crossover m, where the Mott phase is characterized by the gap A m f = /j — J (for 0 < J < p) {26!. 
Going beyond the mean held approximation and using the dimensionless parameters 


cflJ 'ye^ J cflJe^ J 

7 - 1 + e^ “ ~~ 1 — 7 “ 1 + - cflJ 


(5) 


we find the following condition for the existence of a Mott phase: 


Mott correlations: For n > 0 and 0 < A < 1 the exciton correlation function decays exponentially 
as 0 < (e^ tH ale~^ tHH a r ’e~ H ) < Coe~^ r ~ r ^ with a finite prefactor Cq and a finite decay length 
£ < -l/log[(l + e /3J ) 7 ]. 


The derivation of the upper bound is obtained from a hopping expansion, which is given in App. [A] 
It should be noted that this expansion consists of two contributions: In terms of Feynman’s functional 
path integral, the hopping walks of the excitons may or may not cross the time boundaries. The former 
is essential for the derivation of the Mott condition A < 1, whereas the latter requires only the weaker 
condition 7 < 1 . 
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3 Discussion 


The interaction in the exciton gas is mediated by the Pauli principle of the constituting fermions (electrons 
and holes), leading to a hardcore interaction of excitons. This repulsive interaction can stabilize a lattice 
commensurate state, where each lattice site is occupied by one exciton. Without hopping (i.e., for J = 0), 
such a state exists if the the lattice density p = 1/(1 + e -/3/i ) is p = 1. This requires /3p ~ oo; i.e., a 
positive chemical potential at a vanishing temperature. We have used in our discussion that the Mott 
state is characterized by an exponentially decaying exciton correlations. This characterization is weaker 
than to fix n = 1, allowing for gapped thermal exciton excitations. It leads to the condition A < 1 or, 
equivalently, 

c/3J(l + eP J ) < 1 + . ( 6 ) 

This condition reflects the fact that an exciton superfluid state with long-range correlations can always be 
destroyed by strong fluctuations for certain values of the hopping rate (3 J and the chemical potential f3p. 
As an example the phase boundary with A = 1 on the square lattice (i.e. for c = 4) is depicted in Fig. [ 2 j 
Although the destruction of a superfluid phase seems plausible in the case of suppressed tunneling, i.e., 
for /3J -C /3p, the sufficient condition (0 for the appearence of short-range phase correlations requires a 
calculation, for instance, in terms of a hopping expansion (cf. App. 1X1) . 

The temperature enters the condition for an exponential decay only through the normalization of 
the chemical potential and the hopping rate as ftp and /3J, which reflects that either the thermal fluc¬ 
tuations at high temperatures or the quantum fluctuations at low temperatures are responsible for the 
exponentially decaying correlations. As /3J is a measure for tunneling (i.e., quantum fluctuations), in the 
case of /3 J 1 quantum fluctuations dominate. Then we must provide a sufficient large f3p to obtain 
Mott correlations: /3p > j3J. This agrees with the mean field result of Ref. [2(5], although the mean 
field approximation overestimates the stability of the Mott phase against fluctuations for most values 
of p. In the high temperature regime with /3 J -C 1, where thermal fluctuations are dominant, there is 
no restriction for j3p 1 which even can vanish. This is also found in the plot of Fig. [2] Thus, in both 
asymptotic regimes a Mott phase exists. In particular, it is possibles to obtain a Mott phase for p = 0, 
where the density is p « 0.5. This situation should be accessible for gated coupled quantum wells. 

In conclusion, as an extension of our previous consideration in Ref. 20] for a transition to a Mott 
state at high densities, we proved, using a method which goes beyond the mean field approximation, that 
a bosonic Mott phase exists in an electron-hole bilayer through the formation of indirect excitons. In 
this strongly correlated phase, strong fluctuations destroy the long range correlations in the exciton gas 
at any temperature, either by thermal or by quantum fluctuations. 


Appendix 


A Functional integral representation 


To prove the existence of Mott correlations we consider a Grassmann functional integral representation 
of the partition function Z = Tre~^ H for space-time variable 
1, 2,..., M [23]: 


Z= exp ('IpXxfoxfoxi’Xx + ^lx^‘2xVxx'i>2x'i’lx') , Vrn,r 

Jib 


= (r,n) 

with the discrete time n = 

1 + mV 

for r' = r, n! = n + 1 

ij , 

M ° rr 

for r' 7 ^ r, n! = n + 1 

0 

otherwise 


with the Grassmann integration /]. [23] . At the end we take the limit M —> oo. It should be noticed that 
the Grassmann variables ip x are anti-periodic in the time direction. 

The correlation function of Eq. ([3]) reads in terms of the functional integral 

^ly^Xv^Xy'^ly') = 7? J ^yi’Xyi’Xy'^y' exp (lpXx4’2xi>2x4’Xx + ^Xx^2 x V X x' (8) 
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Figure 3: Walk from y = to y' = ( r',m '). The thick line represents the connecting walk, each 

dashed line has the weight 1 + and the dotted line weight 1. There is either a walk that connents y 
and y’ directly (upper example) or through the periodic boundaries in time direction (lower example). 
The weight of the walk is estimated in (THl) . 


with y = ( r,Mt) and y' = ( r',Mt'). It is convenient to introduce the generating functional 

Z {{ a VV'}) = / exp + V , l®V’2x( u + Ot) xx >'ip2x''*Plx') 

Jlf> 

and take the derivative of log (Z({a yy r}). This, for instance, provides for the correlation function 


(9) 


(foylplylply'foy') = 


d 


da u 


log (Z({a yy >}) 


Q !=0 


A.l Mott phase 


The Grassmann integral © can be directly calculated if J rr > = 0 (i.e., in the absence of exciton hopping), 
giving Z M = [1 + (1 + (3n/M) M ] N for a lattice with N sites. For M —> oo we obtain limM^^ Z M = 
(1 + eP^) N . In this case the exciton density is p = 1/(1 + e -/3/i ). This suggests that we consider the 
hopping term 


'^lx' , p2x j^jWxx 


('Wxx' — Jrr'3n' ,n+l) 


( 10 ) 


as a perturbation and apply the linked cluster approach for log Z{{a XX '}), which leads to a connected 
walk with end points y and y' |27j (cf. Fig. [3]). In other words, we expand in (|8|) 


1p2yi>l y tply'i>2y' e n exp(i’lxi’2xV X x>ll’2x'll>lx') 

W ( x,x')EW 
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( 11 ) 


= tp2yi>lvi>ly'i>2y' e n (1 + i>lxi>2xV xx 'i>2x'^\x’) 

W (x,x')EW 

along the walk It 7 between y and y'. Here it should be noticed that W can either connect y and y' directly 
(upper example in Fig. [3]) or by using the Grassmann anti-periodic boundaries in time (lower example 
in Fig. 0 . The latter case will be called a walk with boundary crossing. 

After the Grassmann integration expression m gives for walks with k boundary crossings a sequence 
of m — nn! + Mk factors v xx r. 

{4 > 2y'4 , ly4 , ly'' l l , 2y l ) — ^ ' (1 + Vv yxi v xiX 2 ' ' ' V x m _ m , +Mk y' j (12) 

fc>o x 1 ,...,x m _ m , +Mk ew 

where n v is the number of visited spatial sites. The product of the v terms for a fixed walk W from 
z = ( r,m ) to z’ = ( r',m') with 0 < m' <m < M and without boundary crossing is estimated as 


0 7 Vzxi’^xix 2 * ' * Vx n 


< 


(1 + 


1 + m(M + J) 


m—m'—n v -\-1 


m 


n v —1 


n = 1 


n > 2 


(13) 


Now we use the time variables t = m/M, t' = m'/M and take the limit M oo to obtain 


0 < (1 + e^) nv v ZXl v Xia 


Vx ™-™' z ' ~ \ \/3J/( 1 + e ^]n.-'i e ^+ J )( t - t ')/( 1 + e ^) n v >2 

e /3(^+J)(t-i')/( 1 + e /3#x) < < e /3J ^ 


J e ^(t-i')/(i + e ^) 


n v = 1 


Moreover, we have 


■ (14) 


(15) 


since 0 < t — t' < 1 . As a time independent upper bound we can choose t = 1, t' = 0. 

Finally, we must sum over all possible walks W. There are c (number of nearest neigbors) choices for 
each hop of an exciton. With k — 1 boundary crossings this implies for the correlation function 


{4’2yi’ly^ly'4’2y') < >) «nQn(| 


n> 1 


E E < 

k> 2 ni,...,rifc>2 


7 rifc 0 rai, 


.(k -r'\) 


(16) 


where 


®ni,...,n k (\r - r'\) = | J 


for ni + ••• + rife — fc > \r 
otherwise 


r’ | + 1 


a n 


1 n = 1 

e PJ^n -1 n > 2 


and 7 is defined in Eq. ([5]). & ni ,...,n k (\r — r’\) enforces that at least \r — r'| + 1 sites are visited to connect 
the sites r and r'. Next we perform the summation in Eq. © without the factor Q ni ,...,n k {\i —r'|) and 
obtain 


e«-+e e 




= l + e / 3 J ^ 7 n - 1 + ^e fe/3 ‘ / = 


1 + (e l3J - 1 ) 7 A 2 


i>2 


k> 2 


2 > 2 


n> 1 /c>2 ni,...,nfe>2 

provided that A = 'ye ^ J /(I — 7 ) < 1. This condition can also be written as 


1-7 


c/3J(l + e pJ ) < 1 + 




1 — A 

(17) 

(18) 


Thus, the contribution of the summation of crossings is a factor 1 + e@ J . Finally, we must include the 
factor 0ni,...,n m (l 7 — r'|) in the summation to obtain the upper bound in (ITBl) . This is equivalent of taking 
in m only terms into account with powers n p r r I + 1 and higher. The first term on the right-hand side 
of Eq. (11711 is a geometric series in powers of 7 and the second term a geometric series in powers of 
7(1 + e^ J ), since 


A 

1 - A 


7e /3J 

1 — 7(1 + e@ J ) 


= 7e / 3 J 5 >( 1 + e /3J )]' • 

l> 0 
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Thus, the second term gives the leading contribution in Eq. m, which allows us to extract 


(foylplylply'foy') < C 0 


c/3J( 1 + e^ J ) 
1 + ePv 


\r'—r\ 


as an upper bound for the correlation function (1161) . 


(19) 
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